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ABSTRACT 

High redshift supernovae and Cosmic Microwave Background data are used to con- 
strain the Cardassian expansion model (Freese & Lewis, 2002), a cosmology in which 
a modification to the Friedmann equation gives rise to a flat, matter-dominated Uni- 
verse which is undergoing a phase of accelerated expansion. In particular, the precision 
of the positions of the Doppler peaks in the CMB angular power spectrum provided 
by the Wilkinson Microwave Anisotropy Probe (WMAP) tightly constrains the cos- 
mology. The available parameter space is further constrained by various high redshift 
supernova datasets taken from Tonry et al. (2003), a sample of 230 supernovae collated 
from the literature, in which fits to the distance and extinction have been recomputed 
where possible and a consistent zero-point has been applied. In addition, the Cardas- 
sian model can also be loosely constrained by inferred upper limits on the epoch at 
which the Cardassian term in the modified Friedmann equation begins to dominate 
the expansion (zgg). Using these methods, a Cardassian cosmology is constrained at 
the 2tT level to 0.19< ri„i<0.26, 0.01< n <0.24 for the Cardassian expansion param- 
eter, n, and 0.42< Zf,q <0.89, in contradiction to the previous constraints of Sen & 
Sen (2003a). There is also a large discrepancy between the la confidence regions de- 
fined by the CMB and tightest Supernova constraints, with the CMB data favouring a 
low-f2m, high-n Cardassian cosmology and Zgg >1, as opposed to the Supernova data 
which supports a high-r2„j, low-n cosmology. 



1 INTRODUCTION 



Over the last ten years, the m-z relation below z ~1 has 
been constrained through observations of distant type la 
supernovae, and indicates that the Universe is currently in 
a phase of accelerated expansion (e.g. Perlmutter et al. 1997; 
Riess et al. 1998). This is supported by recent measurements 
of the Cosmic Microwave Background angular power spec- 
trum and the galaxy power spectrum. The standard model 
provides an interpretation of these results by proposing that 
the Universe consists of only ~4 per cent baryonic matter 
with the remainder in the form of Cold Dark Matter and 
a dark energy component which drives the expansion. This 
cosmology provides a good fit to the Wilkinson Microwave 
Anisotropy Probe (WMAP) CMB power spectrum and high 
redshift supernova data, and is consistent with local large- 
scale structure and cluster baryon fraction measurements of 
the matter density. This cosmology is slightly problematic 
however as it invokes the use of two unobserved forms of en- 
ergy. While the evidence for Cold Dark Matter is consider- 
able, evidence for the existence of a dark energy component 
remains indirect. 

An interesting alternative to the standard model of cos- 
mology recently proposed by Freese & Lewis (2002) invokes 
an additional term in the Friedmann equation, proportional 
to p", which results in a Universe which is consistent with 
the recent evidence for an increasing expansion rate, but is 



both flat and matter dominated. The need for a dark energy 
component is therefore removed, and the expansion is driven 
solely by the new p" term in the modified Friedmann equa- 
tion. The theoretical motivation for this Cardassian model 
of cosmology (see Freese & Lewis 2002) is fairly specula- 
tive however, and the implications for Einstein's equations 
of General Relativity are currently undetermined. 

In order to constrain the Cardassian model, it is criti- 
cal that the observational data used do not assume a prior 
cosmology. The simplest constraints can therefore be im- 
posed by the CMB angular power spectrum, high redshift 
supernova data and local large-scale structure, since in a 
Cardassian cosmology these observables depend on Q.^ , the 
Cardassian parameter n, and the redshift at which the Car- 
dassian term in the modified Friedmann equation begins to 
dominate, Zeq. 

Previously, Sen & Sen (2003a) used predictions for the 
locations of the first and third Doppler peaks in the CMB 
power spectrum for a Cardassian cosmology, and Archeops 
and Boomerang data to constrain possible values for Qm 
and n. Supernova data from the Calan-Tololo project and 
the Supernova Cosmology Project (Perlmutter et al. 1997) 
were also used to constrain these parameters, since the mod- 
ification of the luminosity distance in a Cardassian cos- 
mology alters the predicted apparent magnitude for super- 
novae. These two constraints restricted the parameters to 
0.31<n<0.44 and 0.13< 0.23. Zhu & Fujimoto (2003) 
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Figure 1. Contours of Qm with the Cardassian parameters n and 
Zeq- n^=0.1 (dot-dashed), 0^=0.2 (dashed), f2m=0.3 (solid) 
and Qto.=0.4 (dotted) tracks are shown. 



Figure 2. Contours of Zeq with the Cardassian parameter, n, 

and Qm- Contours for Zeq=1.2 (solid), Zeq=l.O (large-dashed), 
2e(j=0.8 (dotted), Zeq=0.6 (small-dashed), Zeq=OA (dot-dashed) 
and Zeq=0.2 (dot-dot-dashed) are shown. 



also used these supernova data to constrain the Cardas- 
sian model, and found best fit parameters of n=-1.33 and 
Zeq=0.43 assuming 0^=0.3. However, it is not clear how the 
fit is affected for different values of flm- More recently Sen 
& Son (2003b) used WMAP and Boomerang measurements 
of the locations for the first, second and third Doppler peaks 
to further constrain the parameter space. 

In this paper, the parameter space is constrained using 
the WMAP and Boomerang measurements for the locations 
of the first, second and third Doppler peaks in the CMB an- 
gular power spectrum used in Sen & Sen (2003b), and vari- 
ous high redshift supernova datasets taken from Tonry et al. 
(2003). An orthogonal constraint is also applied using an as- 
sumed upper limit of Zeq — l; this value is used in order to be 
consistent with local large-scale structure (Freese & Lewis, 
2002). In Section 2, the Cardassian model is outlined. The 
supernova analysis method and the resulting constraints for 
various samples are presented in Section 3, along with the 
CMB constraints and the additional constraint imposed by 
the inferred upper limit on Zeq- The discussion and conclu- 
sions follow in Section 4. 



at Zeq, the epoch at which the Cardassian term begins to 
dominate. Since /?m=Pm,o(l+z)^ 

B = ^[p^,oil + Zeqf]'~" (2) 

assuming the radiation contribution to the energy density is 
negligible. This expression gives rise to a modification in the 
critical energy density such that in a Cardassian cosmology 

Pcrit = F{n, Zeq) X Pcrit,old (3) 

where 

Fin,Zeq) = [l + il + Zeqf^'-''^]-' (4) 

which for a flat cosmology is equivalent to fim ■ Hence, this 
modification to the critical energy density is consistent with 

a flat, matter dominated cosmology. Fig. 1 shows the vari- 
ation of fim (=F(n,-2:eq)) with n and Zeq- Equivalently, Fig. 
2 shows the variation of Zeq with Qm and n. 



2 THE CARDASSIAN MODEL 

An additional term is added to the Friedmann equation such 
that 

H^ = Ap + Bp^ (1) 

where A=^^, H is the Hubble constant, p is the energy 

density and n is a parameter of the Cardassian model. When 

2 

the second term dominates, the scale factor a (x 1 3n , gives 
rise to accelerated expansion for n < |- In order to deter- 
mine B, it is convenient to equate the first and second terms 



3 OBSERVATIONAL CONSTRAINTS 
3.1 Supernova Data 

In order to constrain the cosmological parameters of the 
Cardassian model, firstly, high redshift supernova data is 
used. Since altering the cosmology modifies the luminosity 
distance, the Cardassian model can be constrained by com- 
paring empirical data with predicted values for the appar- 
ent magnitudes of standard candles at known redshifts. The 
relation between apparent and absolute magnitude can be 
written in the form (Perlmutter et al. 1997) 

mB=MB + ^\ogVL (5) 
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Figure 3. Constraints from all 230 supernovae listed in Tonry 
et al. (2003) on Qrn and the Cardassian parameter n. The Icr 
(darlsest shade), 2(7 and 3a confidence regions are shown. 
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Figure 5. Constraints from the 130 supernovae published by the 
High-2 Supernova Search Team (taken from Tonry et al. (2003)) 
on Qm and the Cardassian parameter n. The confidence regions 
are shown as before. 




Figure 4. Constraints from the 54 supernovae published in 
Perlmutter et al. (1997) in Fit C, including supernovae from 
the Calan-Tololo Survey and the Supernova Cosmology Project 
(taJfen from Tonry et al. (2003)), on Urn and the Cardassian pa- 
rameter n. The confidence regions are shown as before. 
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Figure 6. Constraints from the 172 supernovae (for which a red- 
shift limit of z >0.01 and an extinction limit of Ay <0.5 mag. 
were applied where possible) on Qm and the Cardassian param- 
eter n. The confidence regions are shown as before. 



where Mb=Mb — 5logHo + 25 is the Hubble constant- free 
absolute B-band magnitude, and is determined empirically 
for type la supernovae as Ai s=-3.32±0.05 (Perlmutter et al. 
1997) from the intercept on the Hubble diagram. Di=i/orfi 
is the Hubble constant-free luminosity distance. The lumi- 
nosity distance is given by 



dL{z) = 




(6) 



where E{z) in a Cardassian cosmology is given by (Preese 
& Lewis, 2002) 

E\z) = nm{l + zf + {l-^m){l + zf'' (7) 
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Equations 5, 6 and 7 determine a theoretical predic- 
tion for the apparent magnitude (ms(f2m, n)) of a supernova 
at a particular redshift. Hence, Q,m and n are constrained 
through a determination of — Xmm where 



= E 



(ms(nm,n) — msi) 



(8) 



where (Tms is determined by quadrature of the errors on 

log2?L and Mb- For a more robust determination of the 
errors, the uncertainty in Mg and the distance should be 
integrated over the entire parameter space, and an error on 
the redshift could also be included in the determination of 
X^- However, these makes little difference to the resulting 
confidence regions and the determination of the errors de- 
scribed above is adequate for this analysis. 

The observational data is taken from Tonry et al. 
(2003), which presents redshift and distance information for 
230 supernovae compiled from the literature and eight new 
supcrnovac from the Higli-z Supernova Search Team. Since 
the techniques for analysing the observational data vary be- 
tween individual samples of supernovae, the authors have 
attempted to recompute the extinction estimates and the 
distance fitting using the methods of Riess et al. (1998), Jha 
et al. (2003) and Tonry ot al. where possible, as well as apply 
a consistent zero-point for the local distance calibration. A 
complete description of this procedure is given in Tonry et 
al. (2003). 

Four supernova subsamples are selected from the 230 
supernovae listed in Tonry ct al. (2003). Firstly, the entire 
sample of 230 supernovae is used; the resulting Ict, 2a and 
3cr constraints are shown in Fig. 3. Secondly, the 54 super- 
novae included in Fit C in Perlmutter et al. (1997) taken 
from the Calan-Tololo Survey and the Supernova Cosmol- 
ogy Project are selected (see Fig. 4). Thirdly, the 130 super- 
novae originally published by the High-« Supernova Search 
Team are used to constrain the parameter space (see Fig. 5). 
Finally, a redshift cut of z >0.01 and an extinction cut of 
Ay >0.5 mag. is applied where possible, as in Tonry et al. 
(2003), in order to remove objects where an uncertain ef- 
fect is apparent from either local peculiar velocities or host 
galaxy extinction; the constraints arising from the resulting 
172 supernovae are shown in Fig. 6. 
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Figure 7. Constraints on Qrn and the Cardassian parameter, n, 
from the location of the first (solid) and second (dotted) Doppler 
peaks in the WMAP angular power spectrum, and the third 
(dashed) Doppler peak from Boomerang. The 1<T and 2cr contours 
are shown in each case. 
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(10) 



where ao=l, ais=1100 ^ and 

4-3n / 1 ~ f2m,0 



X{a) = I O -h O 



(11) 



For a spectral index of ns=l, the locations of the Doppler 
peaks are related to the acoustic scale by the relation 



Im =lA{m- (t>m) 



(12) 



3.2 The Cosmic Microwave Background 

In order to further constrain the available parameter space 
of the Cardassian model, empirical measurements for the 
locations of the first three Doppler peaks in the CMB an- 
gular power spectrum arc compared with predicted values 
in a Cardassian cosmology. From Doran et al. (2001), the 
acoustic scale can be defined in terms of the conformal time 

, To — Tls /r,\ 

Ia = 7r (9) 

CsTls 

where tq and tu are the conformal time at the present day 
and at last scattering respectively. Cs defines the mean speed 
of sound before last scattering and is constant for a partic- 
ular and is taken to be 0.52 (Doran et al. 2001). The 

Prad ^ ' 

definition derived by Sen & Sen (2003a) for this expression 
of 1a in a Cardassian cosmology is 



where (\>m is a phase shift arising from driving and dissipative 
effects in the photon-baryon fluid before last scattering, and 
is weakly dependent on and n^ (see Doran & Lilley 

2002). The expressions for 0^ for the first three acoustic 
peaks are given in Sen & Sen (2003b). Hence for assumed 
values of Slb/i^ and n^, the locations of the Doppler peaks 
in a Cardassian cosmology can be predicted for a particular 
f2m and n. In the following analysis, values of f26/i^=0.02 
(e.g. Buries et al. 2001; Pettini & Bowen 2001; Kirkman 
et al. 2003) and n^^LO are used, the effects of which are 
investigated in Sen & Sen (2003b). 

The WMAP value for the locations of the first and 
second Doppler peaks are given as li=220.1±0.8 and 
l2=546±10 by Hinshaw et al. (2003) from Gaussian and hy- 
perbolic fits to the power spectrum. The location of the third 
peak is given as l3=825li^ (Hu et al. 2001). The resulting 
1<7 and 2(7 constraints arising from each peak location are 
shown in Fig. 7. 
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Figure 8. Here the constraints from the CMB data, the 172 
supernovae sample (with Ay <0.5 mag. and z >0.01), and the 
upper limit on Zeq are shown. The la, 2a and 3<t confidence 

regions from the first three Dopplcr peaks have been combined 
and are indicated by solid contours. The supernova confidence 
regions arc shown as in Fig. 6. The Zeq=l-0 track is indicated as 
a large-dashed contour as in Fig. 2. 



4 DISCUSSION &: CONCLUSIONS 

The Cardassian expansion model utilises an additional term 
in the Priedmann equation, proportional to p", such that 
the Universe is flat, matter-dominated and in a phase of ac- 
celerated expansion. The need for a dark energy component 
is therefore removed, as the expansion observed in distant 
supernovae is driven solely by this additional component. In 
order to constrain the parameters of the Cardassian model 
n, Qm and z^q (the epoch at which the Cardassian term 
begins to dominate the expansion), it is necessary to select 
observational data free from any cosmological priors. 

In this paper, various samples of type la supernova data 
from Tonry et al. (2003) were used to constrain Q-m and n. 
Fig. 4 shows the la, 2a and 3cr confidence levels derived 
from the 54 supernovae included in Fit C of Perlmutter et 
al. (1997). This is in general agreement with the constraints 
imposed by Sen & Sen (2003a) for the same sample of super- 
novae, although the constraints arc slightly broader, perhaps 
due to their use of contouring for absolute values of . Since 
the errors are uncertain, more robust confidence levels can 
be derived when contouring around a minimum value of 
as used for the supernova data in this paper. Sen & Sen 
(2003a) also use a different method of analysis in which a 
different value of A^b is used with an uncertain associated 
error. 

The tightest constraints from the supernova data are 
provided by the 172 supernovae sample (Fig. 6) for which an 
extinction limit of ^v=0.5 and a redshift cut of z >0.01 have 
been applied in order to remove objects with potentially 
large effects from perculiar velocities at low redshift and 
large host galaxy extinctions. The resulting la constraints 



favour a value of 0.19< Sim <0.32 and n <0.13. However, the 
available parameter space within the Scr confidence region 
is large, although this represents a significant improvement 
over previous constraints. 

To constrain the available parameter space further, 
WMAP (Hinshaw et al. 2003) and Boomerang (Hu et al. 
2001) measurements for the first three Doppler peaks in the 
CMB angular power spectrum are compared to the associ- 
ated Cardassian prediction as in Sen & Sen (2003b) (Fig. 7). 
These three constraints arc in remarkably good agreement 
at fim ~0.15, but begin to diverge at the values of Q.m and 
n suggested by the 172 supernovae sample. 

In Fig. 8, the combined constraints from the first three 
Doppler peaks and the 172 supernovae sample are combined. 
Also included is the Zeq = l-0 track which represents an ap- 
proximate upper limit arising from the effect of the Car- 
dassian term in the modified Friedmann equation on local 
large-scale structure (Freese & Lewis, 2002). The la con- 
straint imposed by the combined CMB data indicates values 
of 0.15<nm<0.18, 0.28< n <0.36 and Zeq >1.0. However, 
there is good agreement at the 2a level with constraints 
of 0.19< Qm<0-26, O.OK n <0.24 and 0.42< Zeq <0.89. 
This is consistent with the loose constraint of Zeq <1.0 im- 
posed by local large-scale structure. However, it is in con- 
tradiction to the constraints derived by Sen & Sen (2003) of 
0.13< 0^^0.23 and 0.3K n <0.44. 

In conclusion, high redshift supernova data and recent 
mesaurements of the CMB power spectrum provide rela- 
tively tight constraints on the available parameter space 
of the Cardassian expansion model. There is a large dis- 
crepancy between the CMB and supernova la confidence 
regions, with the CMB la constraint implying a parame- 
ter space which is potentially at odds with local large-scale 
structure and is rejected by the supernova data at the >2a 
level. However, there is good agreement at the 2a level, with 
the associated confidence regions limiting 0.19< Qm^O-26, 
O.OK n <0.24 and 0.42< Zeq <0.89. 
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